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ABSTRACT 

For typical models of binary statistics, 50-70% of core-collapse supernova (ccSN) progen- 
itors are members of a stellar binary at the time of the explosion. Independent of any con- 
sequences of mass transfer, this has observational consequences that can be used to study the 
binary properties of massive stars. In particular, the secondary companion to the progenitor of 
a Type Ib/c SN is frequently (~ 50%) the more optically luminous star since the high effective 
temperatures of the stripped progenitors make it relatively easy for a lower luminosity, cooler 
secondary to emit more optical light. Secondaries to the lower mass progenitors of Type II 
SN will frequently produce excess blue emission relative to the spectral energy distribution of 
the red primary. Available data constrain the models weakly. Any detected secondaries also 
provide an independent lower bound on the progenitor mass and, for historical SN, show that 
it was not a Type la event. Bright ccSN secondaries have an unambiguous, post-explosion 
observational signature - strong, blue-shifted, relatively broad absorption lines created by the 
developing supernova remnant. These can be used to locate historical SN with bright secon- 
daries, confirm that a source is a secondary, and, potentially, measure abundances of ccSN 
ejecta. Luminous, hot secondaries will reionize the SNR on time scales of 100-1000 yrs that 
are faster than reionization by the reverse shock, creating peculiar HII regions due to the high 
metallicity and velocities of the ejecta. 

Subject headings: stars: evolution - supergiants - supernovae: general 



1. Introduction 



The discovery of the blue progenitor to SN 1987A (iGilmozzi et al.ll 1987b triggered interest in the bi- 
narity of core-collapse supernova (ccSN) progenitors, since mass transfer prov ided a simple explanation for 
how a relatively low mass ~ 1 5— 20Mm star could explode as a blue supergiant ( Podsiadlowski & Josslll989 : 
de Loore & Vanbeverenll 19921 ). Later population studies showed that ~ 25% of masssive stars may undergo 
sufficient mass transfer to transform a primary that would explode as a red supe rgiant in a Type II SN into 
a stripped Helium (o r beyond) star that would explode as a Type Ib/c SN (e.g., IPodsiadlowski et al.ll 19921 : 
Eldridge et al.ll2008r) . If the secondary gains enough mass as part of this process, its evolution accelerates 
and it can explode as a Type II SN while still a blue supergiant. 
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Many aspects of these possibilities were confirmed by the Type lib SN 1993 J in M81 (3.6 Mpc). A bi- 
nary companion w as first suggested bas ed on excess blue emission in the spectral energy distribution (SED) 



observations (Ma und et al.l 2004 



of the progenitor (Alderin g et al. 119941). and later confi rmed in post-explosion imaging and spectroscopic 



Maund & Smarttll2009h. Binary evolution models produce SN 1993J start - 



ing with two similar mass stars ( Podsiadlowski et al. 1993 : Maund et al. 2004 : Stancliffe & Eldrid~gell2009h . 
Mass transfer removes almost all the Hydrogen envelope of the primary prior to the explosion, leading to 
the Type lib transition from a Type II to a Type lb spectral type as the SN evolves. 



While dir ect searches fo r SN progenitors are now common, direct constra ints on binary 



Li et al. 



companions are 



20071 : iGal- Yam et al. 



relatively rare. Smartt (2009) reviews the status of progenitor searches (e.g., 
20071 : ISmartt et al.ll2009l) . For Type IIP SN there are 20 events with adequate pre-SN data, with 6 detections 
and 12 strong upper bounds. These progenitors are red supergiants with 7M ^ <, M p <, 17Mp, wh ich is 
surprising given that locally red supergiants are found up to masses of 25M (|Levesque et all 12005). No 
progenitor of a Type Ib/c SN has been identified out of 10 cases with adequate pre-explosio n data which 



may be further evidence that mass transfer is playing an important role in creating these SN (Smartt 2009). 



Because the data on SN progenitors rarely have broad spectral coverage and there is little adequate 
post-outburst data, there are only a handful of constraints on co mpanion stars othe r than those for SN 1993 J. 
Two additional Type lib SN, SN 2008ax icrockett et all boosh and SN 2001ig (fevder et al.ll2006k show 
evidence for a secondary. For SN 2008ax, there is a blue excess to the SED of the progenitor, while for 
SN 2001ig there is a candidate for the surviving, blue secondary in post explosion im ages. The existence of 
a companion is limited to M s <, 20M & for the Type Ic SN 2002ap flCrockett et al.ll2007l). which also h as the 
tightest limits on the progenitor of any Type Ib/c SN, to L ^ 2L^ for SN 1987 A (IGraves et all 120051) . and 



to roughly My ^ 2 mag for the Type lib SN Cassiopeia A (IFesen et al. 



the Type IIP SN 2008bk (IMattila et al.ll2008l) and SN 2005cs (IMaund et al 



2006 



Krause et al. 



2005 



2008J). Finally, 



Li et al.ll2006r) have upper 



bounds on excesses to their SEDs at V-band of approximately 1 mag. That the constraints are few is not 
surprising, given that it will generally be more challenging to identify the secondaries. We should note 
that there is also indirect e vidence fo r the presence of a binary companion from the structure of some radio 
supernova light curves (see lVan Dykll2004l) . 



Five years ago, limits and detections of SN progenitors would have been equally limited, but the passage 
of time renders all studies of SN progenitors and their secondaries easier, since both the number of nearby 
SN and the amount of useful archival data continuously increase. This suggests that an examination of our 
expectations for the properties of binary companions to ccSN progenitors and of any strategies for detection 
are worthwhile. In this paper we first estimate the fraction of ccSN that have a stellar binary companion 
at the time of the SN in ^2] and then examine the optical properties of these companions in $3] In $4] 
we examine whether these models are constrained by any of the existing data on SN progenitors and their 
companions. In $5]we point out that there is a simple technique to discover and unambiguously confirm the 
brighter secondary companions to ccSN even when the position of the ccSN is poorly constrained. Finally, 
in $6] we summarize the results and discuss the observational prospects. We examine these questions in the 
"wide binary" limit, where there have been no interactions between the two stars. In most cases, adding the 
interactions will enhance the visibility of the secondary, as mass transfer frequently leads to a hotter primary 



-3- 



radiating less in the optical and a more massive and luminous secondary. We will consider these effects in 
Kochanek et al. (2009, in preparation). In an Appendix, we discuss other possibilities for locating historical 
supernovae sufficiently accurately to search for surviving secondaries. 



2. The Fraction of ccSN in Stellar Binaries 

We assume that the distribution of primary masses M p is Salpeter, dN/dM p oc M~ x with x = 2.35, and 
that fraction F of the primaries have binary companions distributed in mass M s as f{q) with q m [ n < q = 
M x /M p < q max where J dqf(q) = 1. With these definitions, the joint distribution of primary and secondary 
stars in mass is 

—^L — = FM-„ x - l f(q) , ( 1 ) 

dM p dM s p JKHh K ' 

and the distribution of secondary masses is 
where 

rqmax 

f q = / tf^mdq. (3) 



mm 



Note that the distribution of secondaries in mass is also Salpeter in form, and that at all masses there are 
Ff q secondaries for every primary. When we observe a ccSN, it can be the collapse of a single star that was 
never in a binary, the primary of a binary or the secondary of a binary, and the relative probabilities of these 
three cases are 

1 - F F F f 

fiineie = and f„ = and = — . (4) 

J 8 l+Ff q Jp l+Ff q l+Ff q 

As long as the primary mass function is a power-law and the secondary mass function depends only on q, 
these fractions are mass-independent. The fraction of ccSN in stellar binaries is simply f p . 



We consider three models for the distribution of secondaries. The first two are based on lKobulnicky & Fryer 
(120071) and are limited to 0.02 < q < 1. The "low mass" model, f{q) oc q~ 6 , is biased to low mass binary 



companions, while the "u niform" model with f(q) cons tant gives equal weight to all masses. The third, 



twins" model is based on IPinsonneault & Stanekl (120061) and puts half of the companions in a twin distri 



bution with 0.9 < q < 1 and the other half with 0.02 < q < 0.9, using uniform distributions for each mass 
ratio interval. The fractions of stars that are secondaries at a fixed mass if F = 1 are f q = 0.288, f q = 0.434, 
and f q = 0.644 for the low mass, uniform and twins models, respectively. Table Q] summarizes the fractions 
(Eqn.HJ) of ccSN that will be isolated stars, binary primaries and binary secondaries for each of these models 
and binary fractions of F = 1, 3/4, 1/2 and 1/4. Biasing the binary distribution to lower masses increases 
f p because there are fewer, diluting massive secondaries. Reducing the binary fraction F, particularly for 
high binary fractions, has a weaker effect on f p , because the change is partly compensated by the reduce d 



numbers of exploding secondaries. If the binary fraction F ^ 0.7, as favored by lKobulnicky & Fryer! (|2007|) . 
then 50-80% of ccSN are part of a stellar binary at the time of the explosion. 
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Table 1. Binary Status at Time of ccSN 



Model 


F 


f single 


f P 




low mass 


1.00 


0.00 


0.78 


0.22 




0.75 


0.21 


0.62 


0.18 




0.50 


0.44 


0.44 


0.13 




0.25 


0.70 


0.23 


0.07 


uniform 


1.00 


0.00 


0.70 


0.30 




0.75 


0.19 


0.57 


0.25 




0.50 


0.41 


0.41 


0.18 




0.25 


0.67 


0.23 


0.10 


twins 


1.00 


0.00 


0.60 


0.40 




0.75 


0.17 


0.50 


0.33 




0.50 


0.38 


0.38 


0.25 




0.25 


0.64 


0.21 


0.14 



Note. — The Model column defines 
the f(q) distribution of the secondary stars 
(see text), F is the fraction of primaries 
that are in binaries and f single, fp and f s 
are the fraction of ccSN that are isolated 
stars, binary primaries and binary secon- 
daries, respectively. 
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Table 2. Observability of Binary Companions 



Model 


M p 


Bband 


I band 


I s < 26 


I s < 26 


I, < 24 


Is < 24 




Mq 


> 10% 


> 100% 


> 10% 


> 100% 


(10 Mpc) 


(3 Mpc) 


(10 Mpc) 


(3 Mpc) 


low mass 


10 


0.32 


0.09 


0.02 


0.00 


0.02 


0.15 


0.02 


0.03 




20 


0.26 


0.05 


0.02 


0.00 


0.09 


0.33 


0.02 


0.14 




30 


0.09 


0.02 


0.07 


0.02 


0.17 


0.42 


0.06 


0.23 




40 


0.63 


0.46 


0.65 


0.47 


0.21 


0.48 


0.06 


0.28 


uniform 


10 


0.45 


0.15 


0.02 


0.00 


0.02 


0.25 


0.02 


0.05 




20 


0.38 


0.09 


0.04 


0.00 


0.16 


0.45 


0.04 


0.22 




30 


0.14 


0.03 


0.12 


0.04 


0.28 


0.54 


0.09 


0.35 




40 


0.65 


0.56 


0.66 


0.57 


0.29 


0.58 


0.09 


0.37 


twins 


10 


0.48 


0.34 


0.10 


0.00 


0.10 


0.38 


0.10 


0.22 




20 


0.45 


0.24 


0.16 


0.00 


0.34 


0.48 


0.16 


0.37 




30 


0.34 


0.12 


0.33 


0.19 


0.40 


0.52 


0.32 


0.44 




40 


0.58 


0.54 


0.59 


0.54 


0.14 


0.55 


0.05 


0.23 



Note. — The model column defines the f{q) distribution of the secondary stars (see text), M p is the primary 
mass, the B band and I band columns give the fraction of secondaries exceeding 10% and 100% of the primary's 
flux in that band. The remaining columns show the fraction of SN with secondaries brighter than either I = 26 
or 24 mag at either 10 Mpc or 3 Mpc. These are calculated for F = 1 but can be corrected to lower values of F 
using Table[TJ 
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3. Observability 



That companions are common do es not imply t hat th ey are observable, so we next make a simple 
model of their visibility. We took the IMarigo et all (120081) solar metallicity isochrones, called the most 
massive star along the isochrone the SN progenitor, and then selected secondaries from the less massive 
stars along the same isoch rone. Massive stars reach the main sequence rapidly compared to their post-main 
sequence life times (e.g., Einnecker & Yorkel l2007|) . so we can ignore small differences in the pre-main 
sequence evolution of the binary stars and simply start from the zero age main sequence (ZAMS). We used 
the estimated HST/ACS HST/ACS B (F435W), V (F555W) and I (F814W) magnitudes. We use the same 
binary distributions as in $2] and assume F = 1. The resulting distributions can be scaled to lower binary 
fractions using Table [T] 



We also examined the Geneva (ILejeune & Schaereri 120011) and Cambri dge STA RS (Eldridge et al. 
2008 ) model sequences, but were able to use them only to spot ch eck the results. Leieune & Schaerern2001 ) 



provide end-phase optical magnitudes only for M <^ 25M Q , while Eldridge et al.l(|2008l) provide histories for 
discrete masses that are t oo coarsely sampled to turn into isochrones. Where we can compare re sults between 



(2008) and Eldridge et al. 



Mari go et all (120081) an d Leieu ne & Schaereri (|200ll) . they are similar. The final phases of the lMarigo et al. 



20081) models are generally similar over the mass range 10M^ < M <, 5OM , 
but with significant differences in the onset of extre me mass loss. In the Eldridge et al.1 (|2008r) models, the 
transition occurs abruptly at M ~ 28M Q , while in the lMarigo et all (120081) models th ere is a steady shift o ver 



the range from 2OM < M < 40M Q . Thus , the M p = 30M ? case we construct from lMarigo et all (120081) re 



sembles the M p = 20M & case, while in the Eldridge et all (I2008h models it would resemble the M p = 4OM 



case. 



Figure Q] shows the distributions of secondaries relative to primaries in mass and magnitude, along 
with the integral distribution of secondaries in mass for the F = 1 models, including the dilution by the 
ccSN of secondaries. To scale these fractions to a different F, simply multiply by the ratio f p (F)/f p (F = 1) 
from Table Q] We show the distributions for primary masses of M p ~ 10, 20, 30 and 40M Q and scale the 
magnitudes to a distance of 10 Mpc. For the lower mass cases, the stars are brightest just when they collapse, 
and only a narrow range of secondary masses can approach the luminosity of the primary. Most secondaries 
are fainter, blue main sequence stars. The secondaries of more massive stars are more visible for three 
reasons. First, for the same mass ratio they are simply more luminous. Second, the luminosity differences 
between main sequence and evolved stars of similar mass are smaller. Fin ally, mass loss lead s to primary 
effective temperatures too high to efficiently radi ate in the optical. Fo r the IMarigo et all (|2008l ) models we 
see this only for the M p = 40M Q case, but in the Eldridge et all ( 2008 ) models the M p = 30M o case would 
show the same properties. 

In Fig. [2] we combine the distributions in secondary mass and luminosity to examine the distribution of 
secondary magnitudes relative to their primaries. The M p =10 and 20M Q primaries undergo core collapse 
as red stars close to the peak optical luminosity of all stars on the isochrone and generally have a blue 
companion. As a result, it is difficult for the companions to significantly perturb the combined spectral 
energy distribution (SED) that would be measured pre-explosion other than by creating a blue excess. For 
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the M p = 40M Q case, the secondaries are frequently brighter than the primary for the reasons discussed 
above. In general, however, both stars are blue, so the overall SED will not be strongly distorted by the 
combination. While many of the secondaries are also massive stars, most are relatively faint and can only 
be found directly in relatively nearby galaxies (< 10 Mpc). 

Table [2] summarizes several of these observational points for F = 1. We first examined the contribution 
of the secondary to the blue flux of the combined SED by computing the fraction of ccSN where a secondary 
contributes more than 10% or 100% of the blue flux of the primary. For lower mass stars, where the primary 
explodes as a red star, the secondary generally distorts the SED and appears as a blue excess. For the 
low mass and uniform secondary distributions, the secondary is rarely (5-15%) as bright in the blue as 
the primary, but frequently (25-50%) a significant (> 10%) perturbation. Spectral distortions are far more 
common in the twins model. Contamination at redder wavelengths for th ese lower mass st ars is minimal 
and will have no effect on the inferences about Type IIP progenitors by ISmartt et all (|2009i) . For the high 
mass M p = 4OM case, the secondary is frequently brighter than the primary at all optical wavelengths 
because both stars are hot enough for their optical colors t o be saturated . The intermediate M p = 30M Q case 
is the least affected by the presence of a secondary in the Marigo et all (|2008r) models because the primary 
is both intrinsically luminous and has retained a significant Hydrogen envelope. Nonetheless, many of the 
secondaries are intrinsically luminous and relatively easy to detect after the explosion fades. 

We can a lso u se these models to roughly estimate the fraction of systems which will interact. From the 
Marigo et all (boosh models we can estimate the initial and maximum radius of the stars. We distribute the 
semi-major axes a uniformly in loga (|Kobulnicky & Fryerll2007T) from a minimum radius a m i„ to 10 5 a„j,„ 
where the minimum radius was s et so that the initial stellar radii were 1/2 the initial Roche lobe radius, 
following iBelczynski et all (120081) . With these assumptions, 25-33% of the primaries will fill their Roche 
lobes at their maximum stellar radii. The fraction depends weakly on the dynamic range chosen for the semi- 
major axes, rising to 30^15% if we reduce it by an order of magnitude, and dropping to 20-30% if we raise it 



by an order of magnitude. These ranges are consistent with more detailed studies (e.g. JPodsiadlowski et al. 



19921 : iKobulnicky & Fryer] 120071 : Eldridge et al.l 120081) . and should mean that adding interactions changes 



the details of our picture but not the broad outline. 



4. Constraints 



Most accountings of t he role of binaries in SN a r e based on the distrib ution of supernova types (e.g., 



Podsiadlowski et al. 



1992: iKobulnicky & Frverl 120071 : Eldridge et all 120081) . Suppose stars from M, 



< 



M < M s trip explode as non-stripped Type II SN (IIP, IIL, Iln), and stars with M > M str j p explod e as stripped 



Type lib and Ib/c SN where the fraction of stripped SN is about f str j p ~ 35% and M m ,„ ~ 8M Q (ISmartt et al. 



2009). In the absence of any binary interactions and assuming a Salpeter mass function, producing the 
stripped SN requires the stripping mass scale to be M str i p = f^ p ~ 2.2M„„„ ~ 18M . This is much 
lower than the mass scale M ~ 25M^ beyond which models for massive stars los e their envelopes prior to 
core collapse (|Lejeune & S chaeren 1200 ll : Eldridge et al. 2008: iMarigo et al.ll2008|) . although it is curiously 
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similar to the upper bound on Type IIP progenitors found by lSmartt et al.l (120091) - If we now allow fraction f\ 
of the progenitors with M < M str j p to interact and become stripped supernovae by mass transfer processes, 
then the interacting fraction must be = 0.08, 0.17 and 0.22 for M vmp = 20, 25 and 30M Q in order to 
produce the observed numbers of stripped SN, and the fractions of stripped SN created by binary interactions 
are 16%, 38% and 52%. While this is overly simplified, it is roughly consistent with the detailed models 
of binary evolution, and, if only ~ 25% of binaries are interacting, it requires the overall binary fraction 
F to be high, consistent with lKobulnicky & Fryerl (|2007l ). Note, however, that creating too large a fraction 
of stripped SN by binary interactions of lower mass stars may by inconsistent with the strong positiona l 
correlations between Type Ib/c SN and Ha emission from young massive stars (lAnderson & Jamesl l2008). 
Stripped ccSN due to binary evolution are likely to show the weaker positional correlations of Type II SN 
with Ha emission. 

The failure to detect Type Ib/c "progenitors" provides our first constraint. Suppose our M p = 40M Q 
model is typical for massive stars stripped by winds to beco me Typ e Ib/c S N, so that the progenitor detec- 
tion limits are also strong limits on secondaries. Based on ISmarttl (120091) 's upper limits on 10 Type Ib/c 
progenitors of Mr = -4.3, -5.7, -5.7, -6.0, -6.0, -6.1, -6.7 -6.9, -7.0 and -7.3 mag, we can estimate 
the expected number of secondary detections for any primary mass and binary model. The constraint is 
relatively weak, but starting to constrain the binary fraction from above. For M p = 40M Q and F = 1 , the 
expected numbers of detectable secondaries in this sample given the magnitude limits are ~ 0.7, 1.0 and 0.5 
for the low mass, uniform and twins distributions respectively, corresponding to 30-50% probabilities of no 
detections. In this case, the twins distribution has fewer detectable secondaries than the uniform distribution 
because many of the twins are also optically faint stripped stars (see Fig.Q]). The results for any of the other 
secondary distributions are similar - only cases like the F = 1 , M p = 30M Q twins distribution can be ruled 
out with high (95%) confidence. This stronger constraint for 30M Q arises from the lack of stripped secon- 
daries at 30Mq compared to 4OM seen in Fig. [TJ Ho wever, a continued fai lure to detect "progenitors" to 
Type Ib/c, particularly with limits as strong as those of ICrockett et all (120071 ) for SN 2002ap, will begin to 
strongly constrain secondary distributions. 

The existence of some binaries does constrain the models from below, although the limits will depend 
much on the definition of the problem. Consider one example. Assu me there have been 22 nearby SN II 
for which adequate pre-supernova imaging data exists - the 20 SN from lsmartt et al. (2009) plus SN 1987 A 
and SN 1993J. In our present model, where we have no binary interactions, both of these latter SN would 
be normal Type II SN. Based on the binary models of SN 1993J, we compute the number of systems out of 
these 22 that would have a secondary within 10% (50%) of the mass of the primary. To have a 50% chance 
of at least one such system requires F > 0.73, (0.11), 0.36 (0.07) and 0.07 (0.05) for the low mass, uniform 
and twins models, respectively. Alternatively, if we interpret SN 1987A as the SN of a secondary, then we 
must have F > 0.11, 0.06 and 0.04, respectively, to have a 50% chance of having observed such an event. 
The constraints are again fairly weak, except that the binary fraction must be very high for the low mass 
binary distribution to have any likelihood of a system like SN 1993 J. 

Finally, we can consider the absence of V-band excesses at the level of approximately 1 mag for the 
Type IIP SN 2008bk and SN 2005cs, both of which have estimated progenitor masses near M p ~ 10M Q . The 
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fractions of ccSN with a secondary exceeding this limit are approximately 5, 9 and 31% for the low mass, 
uniform and twins secondary distributions and F = 1. Thus, while two such non-detections do not constrain 
models, only a few additional, stronger or bluer limits are needed to begin limiting the twins distribution. 



5. Identifying and Verifying ccSN Secondaries 



There are two additional challenges for the problem of identifying ccSN secondaries. First, if searches 
are only feasible for D < 10 Mpc, the statistics will be poor and grow slowly unless the search can be ex- 
tended to historical ccSN where precise astrometric positions are lacking. Second, once there is a candidate 
secondary, we need some means of ruling out chance alignments with an unassociated star. We discuss the 
general problem of locating historical ccSN in the Appendix, while here we focus on the solution to both of 
these problems for bright ccSN secondaries. 

Once the ccSN has faded, the developing SN remnant (SNR) enters the ejecta-dominated phase. Most 
of the ejecta is cool due to adiabatic expansion. At the surface there is a thin layer of shocked CSM and SN 
material separated by a contact discontinuity (see, e.g.. lChevaherlll982l : IChevalier & Franssonll 19941) . Thus, 
searches for young SNR are not very successful because they search for optical (emission line), radio or 
X-ray emission (see Appendix) created by this thin sheath, and this is generally not very luminous unless 
the CSM is very dense. 

Cool, largely neutral gas is best found by absorption, an d this was d r amati cally demonstrated by the 
discovery of the remnant of the Type la SN 1885 in M31 by iFesen et al.1 (I19891) . SN 1885 is now easily 
visible as a ~ 0"8 diameter absorption disk against the bulge of M31 at the wavelengths of high o ptical depth 
absor ption lines like Ca I (4600A), Ca II H&K, Fe I (3021, 3720A), and Fe II (2300-2600A) dFesen et al. 
20071) . The same effects will be seen for ccSN. For example, if we take a representative Type II model, the 
12M S15A model from lwoosley & Weaver! dl995h . theNa, Ca and Fe yields are of order 9 x 10" 4 , 1 x 10" 2 
and 5 x 1O~ 2 M , corresponding surface number densities in the remnant of order 



^ r„ ™ „i« /5000km/sV /lOOyears V i 
S ~ {2,20,40} x 10 16 ( J ( J cm" 2 . 



(5) 



for expansion rate v and elapsed time t. In full, self-similar models of the eject-dominated phase, the 
expansion of the contact discontinuity is slightly slower (radius oc ?~ 09 ) than this simple, ballistic scaling 
(lChevarierlll982r) . If we integrate over the SNR, then the material is spread over the full velocity range, so 
the optical depth at line center is of order 



T = TO 



( 5000km/s ^ 3 ( lOOyears ^ 2 



(6) 



where t ~ 300jc, 1 100* and 1800a; for the Na I (5890, 5895 A), Ca I (4579-4585 A) and Ca II (3934, 3968A) 
lines, and x is the fraction in the appropriate ionization state. The Fe I complex near 3800A will have a 
si milar optical depth. Provided x i> 10~ 3 of the material is in the ground or first excited state, as is true of 
the IChevalier & Franssonl (119941) models, the SNR will be optically thick to these absorption lines. 
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The key question is the evolution of the ionization fraction over the first millennium. For the (probably) 
Type la SN 1006, no Fe I absorpti on is observed, bu t there is m odest Fe II absorpti on corresponding to a 
surface density of order 10 15 cm" 2 jFesen et all 1988b . Models bv lHamilton & Fesenl dl988h of the SN 1006 



remnant found that Fe I is largely photoionized after a few 100 years w hile there is still a significant fraction 



of Fe II. The Fe I absorption in the Tycho SNR (SN 1572) is also weak (Iha ra et al] |2007). so it seems likely 
that the low ionization states only survive for a few hundred years in Type la SNRs. There is no comparable 
information for ccSN remnants. The Crab remnant (SN 1054) shows no strong, broad absorption lines 



(|Sollerman et al.ll2000r ) and we could find no results for Kepler's SN 1604 or Cassiopeia A (~ 300 years). 



For the same physical size, ccSN remnants are likely to be in lower ionization states than Type la rem- 
nants because their higher masses and smaller expansion rates mean they have more material to photoionize 
and higher recombination rates. However, we must also consider the ionizing effects of the secondary (or 
any other central source). While a full photo-ionization calculation including estimates of line emission 
is well beyond the scope of this paper, we did e xamine the photoio nization of a homologously expanding, 



homogeneous envelope of pure H at the level of lOsterbrocl J (U974J) . We computed the time scale at which 



: pnotoion 
k Nl974h . 



an ejecta mass of M e expanding at 5000 km/s at the outer edge of the SNR is reionized given an effective 
ionizing flux from the secondary of Q photons/s, assuming case B recombination and a temperature of 10 4 K 
for the reionized material. For M e = 10M© and log Q = 49 (an 06 star roughly), the reionization front reached 
the surface of the SNR after 150 years. The reionization time then scales linearly with changes in M e and 
Q. Thus, even in the presence of a very hot, luminous secondary, there will be an extended period during 
which the bulk of the SNR is neutral, although a hot secondary will typically reionize the SNR before the 
reverse shock does so. For hotter stars, the reionization of He will be significantly more important than in a 
normal HII region because of its relatively higher abundance, particularly if the progenitor star was stripped 
prior to the explosion. 

The discovery of the SN 1885 SNR in absorption against the bulge of M31 is difficult to apply to more 
distant SNR. First, the remnant must be resolved in order to produce an observable signal. The shell radius 
is approximately 

r, = g =0*035 ( V )(_LJ)(^) (7) 
D V5000km/s/ VIOOyearsy V D J 

where D is the distance. When the remnant is unresolved, the fractional drop in flux relative to an adjacent, 
unobscured wavelength is 

AF ( v \ 2 ( t \V 3M p c V 

— -°- 5 ^-H5000km/sJ VlOOyearsJ ) ' 

where fb ac k is the fraction of the local emission behind the SNR and we have assumed an HST point spread 
function (PSF) radius of r ps f = 0'.'05 > r s . The signal is heavily diluted by the background light encompassed 
by the PSF but not behind the SNR. Once the remnant is old enough, t ^ 100(D/3Mpc)(v /5000km/s) years, 
to be resolved (2r s > FWHM = O'.'l), the fractional drop of AF/F = fb ac k is simply the fraction of back- 
ground light for an optically thick transition. Second, the change in surface brightness between a broad 
off-band (500A) and narrow on-band (50A) observation must be detectable. For filters of equal efficiency, 
observations using the optimal 1:3 exposure time ratio for a 10:1 bandpass ratio and assuming perfect image 
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subtraction (off-band-on-band) to produce a difference image of the absorbed flux, the noise in the differ- 
ence image is four times worse than that of an image through the broad filter with the same total integration 
time. Thus, the signal-to-noise ratio for a resolved (r s > r ps f) absorption disk is 

fS\ fback fS\ ( „ f ( ( Jexp_\ 1/2 ,0-0.4(^-22.5) (g) 

UL - 4.2 UL V rj \r psf ) UosJ (9) 

scaled by the continuum signal-to-noise ratio (S/N) cont for the HST WFC3/F438W filter and an exposure 
time of approximately one orbit. Significant detections against the unresolved emission in the disk of a 
galaxy are challenging - one needs old (r s several times r ps f), foreground {fback — 1) systems in high surface 
brightness (/x# <, 22 mag/arcsec 2 ) regions of the disk. 

Disks, however, are very different from bulges, because mean surface brightness has little meaning for 
nearby galaxies on angular scales where a significant fraction of the average surface brightness is due to 
individually detected stars. For a resolved star behind the remnant, the signal-to-noise ratio is 

-) ^ — (-) ^ 10 f -^M V2 lO- 4(B - 24 \ (10) 
Nj disk *-2 \" J com ^ 2500j 

This is a far easier observation, and the probability of detection is largely set by the probability of finding a 
star behind the remnant. This has enormous implications when searching for secondaries to ccSN or trying 
to use seconda ries to localize ccSN becau se the secondaries are always inside the remnant. Moreover, as 
pointed out by lOzaki & Shigeyamal (120061) . in the context of searches for binary companions to the progen- 
itor of the Type la Tycho SN, their absorption signature is unique. Foreground stars show no absorption, 
background stars show both red and blue-shifted absorption and companions show only blue-shifted absorp- 
tion. This means that searches for bright (J> 26 mag) companions can be carried out in finite observing time 
(<: 10 orbits) for the hundreds of years between when the direct emission from the SN/SNR ceases and when 
the remnant reionizes. 



6. Discussion 

The current picture of massive stars and the creation of ccSN implies that a large fraction (<; 50%) of 
SN progenitors are in stellar binaries at the time of the explosion. In many cases, the secondary stars are 
detectable in either pre-explosion or late time observations. In pre-explosion images, the secondary to a 
Type II SN is most observable as a blue excess in the progenitor SED with 10% excesses being relatively 
common (25-50%) and 100% excesses being relatively rare (10%-25%). Excesses become more common as 
the secondary mass distribution is skewed to higher masses. The secondary to a Type Ib/c SN will frequently 
(~ 50%J be the more optically luminous star, because the stripped primary is radiating primarily in the 
ultraviolet. Direct observations of secondaries are, however, challenging, and are most feasible for SN 
closer than ~ 10 Mpc. 



There are, at present, no systematic attempts to survey for secondaries other than lRyder et al.l (120061) 's 
deliberate search for a secondary star to the Type lib SN 2001ig. The data that are available permit a binary 
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fraction F = 1, alth ough the limits on Type ib/c progenitors (iSmarta 12 009) and on color excesses for two 
Type II progenitors ( Maund et al. 2005 : Li et al, 2006 : Mattila et al. 20081) begin to weakly constrain F ~ 1 
for a "twins" distribution of secondaries where 50% of secondaries have 90-100% the mass of the primary. 
As progenitor statistics improve due to the accumulation of both SN and adequate archival data, the data 
will begin to constrain these models. It is worth, however, considering the types of observations needed to 
constrain the secondary population. 



In pre-explosion images, the challenge is to recognize that the SED is a blend of two stars, as lAldering et al. 



(119941) found for SN 1993J. For the Type Ib/c SN, this may be difficult because both stars will typically be 
blue, but for Type II SN the signature will generally be a blue excess to a red star. To date, most detections of 
SN progenitors are at best in two bands, allow ing the determination of a luminosity and either one intrinsic 
color or an extinction (see lsmartt et al. 2009 ). Limiting the presence of a secondary will generally require 
at least three bands. Ideally, a complete three band Hubble Space Telescope (HS T) survey of the ~ 40 
nearby (<; 10 Mpc) galaxies that dominate the local supernova rate (~ 1/year, see iKochanek et al.1 12008|) 
would provide an archival legacy for finding both progenitors and secondaries into the future as well as for a 
broad range of other astrophysical studies. We should note, however, that the recent transients in NGC 300 
and SN 2008S in NGC 6946, whatever their underlying mechanism, show that mid-IR observations are 
also crucial to unde rstanding progenitor systems because in some cases they are self-obscured by dust (see 
Prieto et al] (12008b) and related references). 



The most dramatic possibility in pre-explosion observations is the identification of an eclipsing binary 
as the progenitor system. Very little is known about the time varia bility of progenitor syste ms, with upper 
bounds ~ 0.2-0.3 mag for both SN 1987A over its last ce ntury (see Plofkin & Claytonll2004 . and references 
therein) and for SN 1993J over a 150 day period in 1984 (Cohen et al.lll995r) . Few local galaxies have been 
the subject of the relatively deep, long duration monitoring projects that would be needed to detect eclipses 
in these (super)giant stars. Discovery of even a single example would have dramatic physical consequences, 
not only by making a direct link between binary evolution and a resulting supernova, but because of the 
enormous amount of extra physical information that can be derived from the binar y period and eclipse 
properties. These will, of course, be rare, but we note that for the parameters of the iMaund et al. 
model of SN 1993 J, the probability of it having been an eclipsing binary exceeds 10%, with eclipse durations 
of roughly two years occuring every 10 years. While challenging, the necessary data can be obtained by 
ground based monitoring of nearby galaxi es, as illustra t ed by the eclipsing binary in Holmberg IX, a tidal 
dwarf companion to M81, discovered by Prieto et al. ( 2008 a f) using the Large Binocular Telescope. In 
addition to monitoring any SN progenitors, such data would also inventory all luminous eclipsing binaries 
and variable stars (e.g., Cepheids) as well as providing t he basis for determini ng whether any massive 
stars die as failed supernovae without a dramatic explosion (IKochanek et al.ll2008l) . A single epoch, higher 
resolution HST survey is an important complement to any ground based monitoring. While difference 
imaging methods make it relatively easy to obtain light curves of variable sources from the ground even at 
10 Mpc (the AC signal), stellar crowding means that HST is needed to provide absolute calibrations (the DC 
signal). 



Once the SN has faded, it is relatively easy to find and confirm bright secondary stars of ccSN in 
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galaxies closer than D < 10 Mpc. During its coasting phase, the SNR w ill have very high abso rption optical 
depths in lines of Ca I, Ca II, Na I, Fe I, and Fe II, as lFesen et all (119891) and lFesen et al.1 (12007 ) have already 
used to identify the SNR of the Type la SN 1885 in M3 1, unless something photoionizes the gas. The method 
used for SN 1885, searching for a disk of absorption against the emission from unresolved background stars, 
is possible only when the surface brightness is relatively high (;> 22 mag/arcsec 2 ) and after ~ 10 2 years when 
HST can resolve the SNR. This approach may work b est for identifying h istorical ccSN associated with star 
clusters (e.g. like SN 2004am and SN 2004dj in the ISmartt et all (|2009l) sample). On the other hand, any 
resolved star brighter than <> 26 mag that is either behind or inside the SNR should produce a detectable 
absorption signal. The two cases can be distinguished because background stars will show both red and blue - 
shifted absorption while secondaries will show only blue-shifted absorption (lOzaki & Shigeyamal 12006). 
Moreover, if the search is carried out shortly after the SN fades, contamination by a background star is very 
unlikely given the compactness (milliarcseconds after 10 years) of the SNR. Fortunately, ccSN secondaries 
are generally blue stars at the wavelengths likely to show strong absorption features, although some care 
will be needed for stars with significant B aimer absorption in their atmospheres. 

Existing studies of reionizatio n of SNR by e i ther emission from shocked material near the surface 
or ambient UV b ackgrounds (e.g. IChevalierl 1 19821 : Hamilton & Fesenl 1 19881 : IChevalier & Franssonl 1 19941 : 



Fesen etal 



20071 ) suggest that the bulk of an SNR should remain neutral for several centuries but then 



reionize on time scales <, 10 3 years, although the late-time analyses of this problem have all focused on 
Type la remnants. These time scales seem to be broadly consistent with the limited available data on ab- 
sorption in SNR. The effect of secondary stars on the reionization of remnants seems never to have been 
considered. Crude models suggest that hot, luminous O star secondaries can photoionize an SNR from the 
inside on time scales of 10 2 -10 3 years. This will create a very peculiar cross between an HII region with the 
secondary star as the photoionizing source, a planetary nebula, because of the very high He and metal abun- 
dances, and an SNR because of the high ejecta velocities. Full, dynamical, radiative transfer simulations will 
be needed to c haracterize the observa bility of this phase. The obvious trial scenario, the late time spectra of 
SN 1993J (e.g. lMaund & Smarttll2009L is probably little affected by this process because the secondary is a 
cooler B star. 

The first step in a campaign to identify secondaries is to take all nearby ccSN with accurate positions, 
particularly from post-explosion HST imaging, and reexamine the sites. Ideally, one wou l d start with sys- 
tems for which late time observations already exist (e.g. the 11 ccSN examined by lLi et al.1 (|2002h . although 
these are generally more distant than optimal) to ensure that any source is unlikely to be originating from 
the SN/SNR. The next step is to use the absorption method to try to confirm some of the candidates. Here 
the obvious first cases to try are the candidate secondary to SN 2001ig identified by iRyder et al.l (12006) and, 
once the SN has faded enough to examine absorption in the secondary without contamination from emission 
by the SNR, the secondary star of SN 1993J and (possibly) SN 2008ax. The numbers and properties of the 
secondary population will constrain the role of binaries in the evolution of massive stars and the production 
of the various classes of ccSN. They also provide an independent check of the models used to infer the 
properties of progenitors. In a few cases where the secondary is sufficiently bright, it may be possible to 
examine abundances in the ccSN ejecta. In the early phases, even elements with ejecta masses of ~ 10~ 6 M Q 
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(e.g., Barium) have significant absorption optical depths, although detection may be limited by Doppler 
smearing with stronger lines and inferences may be limited by ejecta clumping. Identifying the secondary 
to a historical SN will also establish that the SN was not a Type la and set a lower bound to the mass of the 
primary. 

I would like to thank J. Beacom, R. Chevalier, A. Hamilton, J. Johnson, M. Pinsonneault, R. Pogge, J. 
Prieto, S. Smartt, K. Stanek, and T. Thompson for comments and discussion. This research made use of the 
Sternberg Astronomical Institute supernova catalogs and the NASA/IPAC Extragalactic Database (NED), 
which is operated by JPL/Caltech, under contract with NASA. 



A. Other Approaches to Identifying Historical Supernovae 



In we introduced a simple means of identifying bright secondaries to ccSN even if the ccSN loca- 
tion is poorly known. This method should work for several hundred years after the SN unless something 
photoionizes the remnant. The absorption method will also work if the SNR lies in front of a sufficiently 
bright background star or if it lies in a region of high surface brightness like a star cluste r. Here we discuss 
other possibilities for identifying the positions of historical ccSN with high precision (|Smartt et all (|2009) 
argue for a goal of orde r 10 m as, albeit at higher typica l distances than 10 Mpc). To provide a sense of 
the problem, iBarth et all (11996) and Van Dyk et al.1 (|1999l ) examined the environments of past SN with HST 
imaging data, but could only localize the SN with positional uncertainties of 1-10 arcsec that made anything 
beyond environmental studies impossible. 

The goal of most searches for SNR in external galaxies has been to carry out more uniform survey s 



than are possible in the Galaxy. Sea rches have used optical emi ssion lines (e.g..lMatonick & Fesenll 19971) . 



radio emission (e.g.. lLacey et al.lll997h . and X-ray emission (e.g.. lPannuti et al.ll2007l) 7The optical and radio 



surveys do not find counterparts to historical supernovae unless they had strong CSM interactions such as 
the examples we discuss below. Specific X-ray searches have been more successful , with X-ray co unterparts 
having been identified for most post-1970 SN0 with the record being SN 1941C dSoria & Perna 2008). In 



gener al, there is little overlap between the sources found by the three approaches (see, e.g., iPannuti et al. 
2007). A further problem is that radio and X-ray detections may not lead to good enough optical astrometry 
in these crowded fields due to problems in matching the reference frames with high accuracy. 

Two additional possibilities are late time light curves and dust echos. The S N themselves fade r elatively 



Li et al. 



rapidly. For example, SN 1987A was fainter than its progenitor by ~ 800 days (|Suntzeff et al Jl991h . 
(|2002h observed 11 ccSN (mostly Type Iln) at late times. The four observed later than 1000 days had 
detections interpreted to be the SN but spanning an interesting range for secondaries (My ~ -4 to -13 mag, 
see Fig. |2]). Extrapolating the light curves using the mean decline rates suggests that in most cases the SN are 
fainter than My ~ -4 mag after 2000 days, On these long time scales, the only internal source of energy is 



'See http://lheawww.gsfc.nasa.gov/users/immler/supemovae_list.html 
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44 Ti which produces a maximum luminosity of L44 ~ 10 3 (M44 / 10 4 M Q ) exp(-f /86years)Lfn i f all the decay 
energy is absorbed and 25% of this if only the positrons are absorbed (IWoosley et al.lll989i) . Combined 



with a slope of only 0.013 mag/year, this will be challenging to detect decades later even if the emission 
is concentrated in the optical/near-IR. Some older ccSN are, however, trivial to find deca des after the event 



becau se of CSM i nteractions. In p articular, br ight radio supernovae lik e SN 1979C (Milisavlievic et al. 



2009k SN 1980K dFesen et al.lll999h . SN 1986J JMilisavljevic et al.ll2008h or SN 1993J Jlvlaund et al.ll2004h 
are still being observed. The general distribution of supernovae from the nearly invisible 44 Ti regime to 
the strong CSM interaction regime is not known. One result of searching for secondaries to all ccSN with 
accurate post-explosion images, would be to clarify these issues and provide strategies for the ccSN without 
such data. 

The last possibility we consid er is searchin g for dust echos. Extraga lactic dust echos h ave been im- 
aged for th e Type Ilpec SN 1987 A dCrottslll988h. the Type la SN 1991T Jsnarks et al.lll999h. th e Type lib 



SN 1993J (lSugerman & Crotts 



2002 
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Type la SN 1995E (lOuinn et all 1200(1 and the Type la SN 2008X (IWang et al.ll2008h . They were gener 



2003b . the Type la S N 1998bu Jcappellaro et aD boOlh. the 



ally found using difference imaging techniques, taking advantage of the frequentl y superluminal p attern 
speeds of the echos. While these detections are all relatively close in time to the SN, Rest et al.l (120051) have 
demonstrated that echos of 400-600 year old supernovae can be seen in the LMC, some of which also have 
superluminal apparent motions (up to v e = 3c). The LMC echos are resolved (A ~ 2"5 across) and can have 
relatively high surface brightnesses (up t o 22 mag/ar c sec 2 ). There have also been ground based attempts to 
find dust echos in more distant galaxies. iBoffi et al.l dl999i) examined the sites of histori cal SN for unusu 



ally b lue emission due to scattered light, finding 16 candidates among 64 nearby SN, and iRomaniello et al. 
(2005) searched for polarized emission near 4 SN in M 83 without success. 



The problem with any search for dust echos is recognizing them against the crowded stellar background. 
The successful detections have all taken advantage of the apparent motions of echos a nd difference ima ging 



Rest et al. 



(2005J) are 



multiple epochs to solve this problem. Here we note that observations comparable to 
possible for older SN in nearby galaxies using HST, particularly if there is dust in the foreground of the SN. 
The outer angular scale of the echos is determined by the elapsed time and the foreground dust distribution, 
with 



r out = 0'.'54 



lOyears 



1/2 



Zdust 



10pc 



1/2 ^3Mpc^ 



1 + - 



ct 



^■Zdust 



1/2 



(Al) 



where Zdust is the distance of the dust along the line of sight from the SN. Thus, the angular scales for nearby 
(D < 10 Mpc) somewhat old (decades t o centuries ) ccSN dust echo s are w ell-suited to HST. Theoretical 
models of echos have been explored by Isugerman ( 2003 ) and Patat ( 2005 ). but we will focus on scaling 
from the observed properties of the many century old LMC echos. The surface brightnesses of the brighter 
LMC echos can be very high, but HST will somewhat under resolve them in one dimension because their 
widths (A ~ 0"04(3Mpc/D)) will be somewhat narrower than HST's PSF. Failure to resolve the echo in 
one dimension reduces its effective surface brightness as 2.5 log A/FWHM where FWHM ~ 0"1 is the 
resolution of HST. Even so, high surface brightness echos at D < 10 Mpc should be relatively easy to detect 
if the stellar backgrounds can be suppressed. The existing HST studies have all used HST pre-explosion 
data to subtract the stellar background and render the echos easily visible. This is not feasible for historical 
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SN, but, the strategy used for the LMC detections will work. The proper motion of an echo with effective 
velocity v e is 



fjL = 0'.'02l( — 

c 



year 



(A2) 



and it must move one resolution element, fiAt = FWHM for difference imaging to be used. The required 
time scale is relatively short 

" ^ (sM^) yCarS ' 



At ~5 



(A3) 



even in the absence superluminal apparent velocities. This can b e tested fo r many SN as part of any lat e 
time search for secondary stars or by revisiting the SN studied by lBarth et all (ll996l)JVan Dyk et al] (11999). 



Li et all ((2002). 
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Fig. 1. — Secondary magnitudes at the time of primary core collapse as a function of the initial mass ratio 
M s /M p . The magnitudes of the primary are the values for M s /M p = 1. The B (blue), V (green) and I (red) 
magnitudes are scaled to a distance of 10 Mpc. The dashed lines show the integral distributions of the 
secondaries in the mass ratio including the dilution by collapsing secondaries for binary fraction F = 1 (right 
scale). From top to bottom at high mass fractions they are the twins, uniform and low mass distributions. 
To rescale these to a different binary fraction, multiply by the ratio f p (F)/f p (F = 1) from Table [TJ For the 
Eldridge et all (120081) evolution models, the 3OM case would resemble the 4OM case. 
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Fig. 2. — Integral distributions of companion magnitudes in the B (blue), V (green) and I (red) at a distance 
of 10 Mpc. The low mass (dotted), uniform (solid) and twins (dashed) distributions are shown for a binary 
fraction F = 1 and the diluting effects of collapsing secondaries are included. To rescale these to a different 
binary fraction, multiply by the ratio f p (F)/f p (F = 1) from TableQ] The vertical lines mark the magnitudes 
of the primary at the time of collapse. Note how the intrinsica lly most luminous M„ = 4OM progenitor is 
optically fainter than the three lower mass progenitors. For the Eldridge et al.1 (120081) evolution models, the 
3OM case would resemble the 4OM case. 



